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Abstract: This paper presents some of the results of analyses of body temperature data of two 


. Species of mammals. The animals were subjected to various environmental conditions with re- 
* spect to variations in periodic light cycle with a period close to 24 hr. New statistical techniques 


are used to study the animals as a function of time and change in environment. These techni- 
ques confirm some previously observed behavior and indicate new aspects of the dynamics of 
adjustment to change in photoperiod. 


L. Introduction 


This paper reports very briefly on the analysis of a large body of data which were 
recorded at Princeton University during the last several years. The data consist of 
temperature and activity recordings, taken at either two- or five-min intervals, on 
several species of small mammals. The animals were isolated in cages which were 
subjected to a variety of controlled photoperiods. 

Tables 1-4 show the total lengths of the data recordings and the various condi- 
tions which were imposed. 

Due to the limitation on the length of this paper, and also due to the fact that a 
full report is currently in preparation, we will not give any discussion of the impor- 
tant questions of the data recording techniques, or the data preparation and com- 
puter preprocessing which were carried out. We will simply state here that it appears 
that the resulting data which were finally analyzed by the techniques discussed in 
this paper were of high quality. 

In section 2 we will give a relatively extended discussion of the statistical techni- 
ques which we used in the study of these data. Some of these techniques have not 
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previously been applied in biology, and, in fact, have only been developed within 
the last few years. 

In section 3, we will give a brief summary of the main results of our analyses. 
This summary is condensed and preliminary in nature. Again, a full discussion and 
further analysis are in preparation. 


2. Methods 


The purpose of this section is to describe the statistical techniques which we 
have used in the analysis of these time series data. Since much of the interest in the 
data concerns evidence of changes in the properties of the mechanism which 
generated the data, some novel statistical techniques were required. In particular, 
we have used techniques which do not require that the process remain invariant 
under time translations. With these techniques we have been able to study the tran- 
sitions which occur when an animal’s environment is suddenly changed. We have 
also, of course, applied the somewhat more conventional techniques of time series 
analysis to appropriately selected portions of our data. In general, we have used 
the techniques which do not depend on long time averages to investigate the appro- 
priateness of forming such long averages. These averages are required, as will be 
discussed in greater detail below, in order to form precise estimates of any indica- 
tion of periodic variation in the data. 

We will first give a brief discussion of the computation of periodogram and spec- 
trum estimates and their uses and limitations when applied to the study of nearly 
periodic phenomena. In the second part of this section we will discuss the compu- 
tation and applications of frequency decomposition by use of complex-demodul- 
ates. The materia! in this later part will probably be less familiar to many. readers. 
An effort will be made to give a fairly self-contained exposition of these data- 
analysis procedures. 


2.1. Periodograms and spectra 

In this section we will discuss the computation of periodograms and spectra 
from the standpoint of the statistical data analysis implications of various possible 
procedures. For other aspects of spectrum estimation, other time series proce- 
dures and relevant computational procedures the reader is referred to the short 
bibliography at the end of this paper. 

First, we should provide some motivation for the use of frequency analysis. Any 
series of data points which may represent the values of some variable and which 


are equally spaced in time may be exactly described by a sum of periodic functions. 


This description in terms of periodic functions may be thought of as a description 
which depends on the set of frequencies, Thus, every set of values of a variable at 
equally spaced time points may be uniquely represented by a set of trigonometric 
functions with specified frequencies and amplitudes. Using the representation in 
terms of frequencies has several advantages: 
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1) Any tendency in the original data for periodic variations to oceur in particu- 
lar frequency ranges may be made more obvious in this representation. 

2) Under appropriate conditions, the coefficients which form the representation 
will be independent of each other. This makes the representation much “simpler” 
in that we may consider the coefficients separately. 

3) Under appropriate statistical assumptions this representation provides a means 
of forming statistically reasonable estimates of the parameters of the process which 
may be generating thé observed data. , 

The formai statement of the above discussion is as follows: define the finite 
Fourier series representation of a set of values of a variable as: 


n 
h, (ws) = Dy exp[-iwt]x,, —1<w<n, (1) 
t=1 


where x, is the value of the variable at time point t, for which we have values for 
t= 1,2,...,2, w is the frequency which may be defined at any point on the’ range 
—a to m, and the complex exponential exp[—i wt] is the. complex-valued periodic 
function which may also be written exp[—iw#] = cos wt — isin wr. exp[—iw#] isa 
periodic function of t with period 2*/<» as can be seen by substitution and use of 
standard trigonometric identities. We can also define the inverse Fourier series 
which gives us the values of x, from h,,(w): 


a—] 
1 
x, > ey. exp [ic,t] h,{w,) : t=1,2,...,0 (2) 
j=0 
where 
a = jafn . 


We can see that these expressions give mutually reciprocal transformations by sub- 
stituting eq. (1) to (2) to yield 


n-1 


1 : Lee 
x,=> p exp [ict] (exp [-ico; ]xq + exp [—ic2;2]x> +... 
T= 


+ exp[—icat]x,+... + exp[iw;n]x,,) - (3) 
For each value of ¢ only the term 
n—1 


zu exp [ico,t] exp[—ico; t]x, (4) 
{= 
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will be non-zero. It will be 
n—l 


Da exp [ica;t — iw, t}]x,=(")x; . (5) 
i=0 


While we may define 4,,(w) for any value of w, it is completely determined — given 
- n data points — by its values at the frequency points given by 


ta, = jan , j= 0,...,n is (6) 


This permits us to decompose a series into at most n+ 1 uniquely defined frequency 
components. The periodogram may be defined as: 


PO) == Winco) : (D 


From eq. (1) it would appear that we could evaluate the periodogram at any value 
of w. However, from the standpoint of statistical interpretation it is very important 
to keep in mind that the periodogram is uniquely defined by its values at the n+ 1 
frequency points given by eq. (6). In fact, a more usual definition of the periodo- 
gram is: 


n 2 
PL; =z 2 x,exp[—2rijt/n]| , 7=0,1,...,2. (8) 
f= 


This definition makes the set of uniquely defining frequencies explicit. It must be 
emphasized that since these frequencies are spaced at intervals of width 1/n the reso- 
lution of periodograms (and spectra) is limited to at very best frequency intervals 
of width not less than 1/n. The implications of this fact will be discussed in more 
detail below. 

Current computational techniques imply that the most direct way to compute a 
sample spectrum is to-compute the periodogram ordinates using a fast Fourier 
transform algorithm and then smooth the periodogram to obtain a consistent esti- 
mate of the spectrum. 

In order to appreciate the statistical consequences of computing a spectrum 
estimate from a sample of data we must consider the following. A sample of data 
can be considered as if there were an infinitely long sample multiplied by the func- 
tion: 


wr=l, 1<r<7r, 


@) 


w,=0, otherwise, 
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where-T is the length of the time average used. If we form a spectrum estimate then 
ci 
fa)=f Lo-P2 Za, (10) 
° 


where L(A) is the Fourier transform of eq. (9), £0) = IK,QUL is the spectrum of 
the infinitely long series, and f,(w) is the sample spectrum. It is important to re- 
member that we must form estimates of f,(w) and f,(w). The quality of these esti- 
mates depends on, among other things, L(A) and K,(A). Using this framework we 
can see that there are ways to improve the quality of the estimates by effectively 
changing L(A). In particular, instead of using the rectangular data window given by 
eq. (10) one can use a window which tapers more gradually to zero at each end. 
This has the effect of producing an L(A) which is smoother and permits much less 
contamination of the estimate at one frequency by variation from another neigh- 
boring frequency. 

Now, if we consider the simple case in which f, (ce) is given by 


f,(w) = K + h5(w,) , (11) 


where K = constant, §(co,) = delta function at (,, and # = scale parameter, we can 
derive the shape of f,(«2) and then discuss the precision with which we can deter- 
mine wy. 
Initially we may take L(A) to be the Fourier transform of w(t). Thus 
sind g7 2 
je (12) 
SAT 


In order to form an estimate of f,(w) using eq. (12) we would have to compute 
the periodogram of the data and use the squares of the periodogram ordinates. It is 
well known that this is a risky thing to do.in general but if # > K (ie., our process 
is known to contain a very large deterministic periodic variation at the unknown 
frequency w,) this procedure may be reasonable. It is important to consider this 
procedure because it results in the highest possible resolution of frequencies-at the 
cost of undesirable statistical properties for the estimates. In this case F,(w) takes 
the form 


ILO)? = (7 


: = ig 
sin(w—a,) 5 ) as) 


(w—w,)} T 
This results in a spectram with the shape shown in fig. 1. 
A variety of measures of the width of a peak such as this have been proposed. As 


reasonable as any, for our purposes, is the half power point given by the distance 
along the x-axis between the two points for which: 


(acest a 


(w—w,)3T 


feo) =Kh ( 


ed a 


(14) 
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Fig. 1. | 
This distance is 
Dee 
2Aw ~ 5, sin 1G), (15) 
where 
Aa = |w—w| 
or 
2a . 
2Aw = 3T° (16) 


As an example, if we consider the data required so that Aw = 1/720 cycles/ day 
(this is a shift of 2 min/cycle for a cycle with period one day), we find that 


T, ~ 750 days . (17) 


Given effectively error-free data and a very low variance at frequencies around 
the supposed delta function we could reasonably say that, given these data, the 
delta function appears to be located at . 


w, t+ Aw. (18) 
However, in this procedure we have assumed the existence of a delta function at 
some constant frequency. We should also consider the possibility that either the | 
delta function changes its position or that the spectrum has a sharp but continuous 
peak. If we allow the possibility of some (random) variation in the (now not quite ; j 


precisely) periodic term in the data then we must use up some of our data in order 
to average over this variability. 

If we assume that the spectrum is continuous and locally smooth then, in the 
procedure which we have just outlined each estimate of the spectrum has aX? dis- 
tribution with two degrees of freedom. If we modify the procedure in order to in- 


BIOLOGICAL TIME SERIES 225 


crease the degrees of freedom the data required for a given resolution increase pro- 
portionally with the degrees of freedom. In our example, if we want estimates with, 
say, eight degrees of freedom and still want to maintain Aw = 2/1440 then we must 
have 


2.2. A frequency resolution nomogram 
. As the length or record required (T) is a linear function of the two variables, de- 
grees of freedom (7) and resolution (Aw), we can present the relationships between 


-. these three variables using a tri-linear graph. 


Fig. 2 presents the relationship 
: n 
T= 750 —-——_. (20 
(Aw/w?) ) 


The axes have the following meanings: 


Label Variable Range Units 


A-a Awfus? -5—128 min per cycle 
B-b n 2—1024 - 
C-c T 23-6000 days 
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Awo/u? is derived from: 


2 2(w2—~) 27 Aw 
2 w? 


Ap =P\—P2 = anf rel 


+ 


(e3) wy F ww), 


and is used since it seems natural to measure resolution in time differences per 
cycle. 

The intersections of the three perpendiculars to the three axes give the values of 
the three variables. Specification of the values of any two variables determines the 
value of the third. 


2.3. Complex-demodulation 

The estimation of spectra as described in the previous section requires forming 
averages over the entire length of the data. This procedure has an interpretation in 
terms of a population spectrum only if it is assumed that the process is wide-sense 
stationary. This may be a reasonable assumption for our data during periods when 
the environment is held constant. It obviously is not reasonable for the investiga- 
tion of behavior at, and for some period after, a change in the environment. In 
order to investigate changes in behavior we have employed complex-demoduiation. 

The basic computational procedure is quite simple. We will state this procedure 
first, then indicate its theoretical motivation, and then show some practical proper- 
ties. 

A complex-demodulate of a time-series, x,, at frequency «, is defined by 


x,(w) = L£(expfiwt]x,), (21) 


where x, is the observed time-series at f= 1,....2, and L( ) is an appropriate low- 
pass filter. If the process which generated x, has a spectrum /(w), then 


f(w) = E(x (w) I?) . (22) 


Therefore it seems natural to define a possible estimate of the sample spectrum by 
n 
1 
fo=— Dr ix orl. (23) 
t=i 


This equation clearly shows the relationship between a spectrum estimate and 
complex-demodulates. The average value of the amplitude of the complex-demodul- 
ate is a spectrum estimate at the frequency w. 

A complex demodulate at frequency <2 may be thought of as a time series 
which contains the variation in the original time series which took place near the 
frequency «. To see how this is true consider a time series: 
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x,=Acoswot, (24) 


define £(c2) to be the Fourier transform of the function £( ) which appears in eq. 
(21). For expositional simplicity we will take 


L(w)=1, lwi<Aw, 
(25) 
L(w)=0, lwl> Aw. 
Then, substituting eq. (24) into (21), we have 
x,(@) = LES exp[icot] A (exp[ico, t]. + exp[—iw, z])} 
= 117A (exp[i(wy + o)t] + expfi(w—w,))} . (26) 
Now if w—w, < Aw and w, +w > Aw using eqs. (25) we have 
x,(e) = +A expfi(w—w,)t] . , (27) 
We can now rewrite eq. (27) in terms of amplitude and phase as 
Ix,(w)I? =34? , 
Im{x,(w)] ; 
= ~1 poeta henley = = 
p(w) = tan Re[x,(a)] {w—-W,)t. (28) 


3. Conclusions 


_ In this section we will give a brief summary of our results and present figures 
which indicate some of the main findings. 
“Our results may be summarized as follows: 

1) The computation of spectra and complex-demoduiates yields results conform- 
ing to previous generalizations, that in nocturnal mammals the circadian period in 
LL is greater than the period in DD (see figs. 3,4,7, 10). 

2) Previous evidence concerning the fact that there is a considerable period of ad- 
justment when the photoperiod is changed is also supported by cur computations. 
However, our evidence suggests that this adjustment period, and new equilibrium 
if any, are substantially more complex than had previously been indicated (see figs. 
8,9, 11). 

3) The complex-demodulates indicate (see figs. 8,9, 11) that there is substantial 
phase shift between the fundamental circadian frequency and its first few harmon- 
ics when a new photoperiod is initiated. In fact, in some cases one of the frequen- 
cies in the set of harmonics seems to lock onto the new photoperiod while others 
fail to do so. This appears to be a very interesting new finding, particularly in its re- 
lationship to the various coupled-oscillator hypotheses which have been suggested. 
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Fig. 11. 


4) There appears to be a persistent overshoot in the adjustment to a new photo- 
period if the new photoperiod is close to 24 hr. This overshoot is typically quite 
small, but it appears in many cases and over quite long time periods. 

In conclusion we would like to remark that the body of data which we have 
been studying appears to be the most extensive, in terms of lengths of records and 
replication within experiments, and the most carefully controlled and instrumented 
that have yet been recorded. Our analyses have already provided new empirical evi- 
dence concerning hypotheses which attempt to explain the sources of circadian 
thythms (particularly the coupled-oscillator hypotheses). We hope that as we and 
others continue to study these data substantial further insights will be obtained. 
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